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ABSTRACT 


In  addition  to  the  pole  series  in  the  singularity  expansion  of  a 
scattered s  response  a  ‘possible  entire  function?*  contribution  has  been 
customarily  included  although  considerable  uncertainty  surrounds  the 
role  of  this  response  component.  Entire  functions  are  frequently  found 
in  the  responses  of  circuits  where  they  represent  forced  components. 

The  results  of  this  investigation  indicate  that  in  certain  extreme 
system  configurations,  like  circuits,  the  forced  response  to  a  time- 
limited  excitation  is  an  entire  function.  In  certain  less  extreme 
configurations,  like  most  scatterers,  the  forced  response  takes  the 
form  of  a  ^quasi-entire  function1*  which  is  composed  of  highly-damped 
pole  contributions.  We  demonstate  that  true  entire  functions  themselves 
can  be  represented  by  infinite  pole  series  whose  poles  have  become 
infinitely-damped.  For  many  scatterers  two  groups  of  poles  can  be 
distinguished:  the  dominant  poles  and  the  highly-damped  poles  whose 
contributions  comprise  a  quasi-entire  function.  The  numerically-derived 
quasi-entire  function  conponents  of  a  conducting  sphere-capped  cylinder 
and  a  conducting  loop  are  determined  by  subtracting  the  dominant  pole 
series  from  the  total  responses  of  these  scatterers.  These  results 
indicate  that  the  forced  response  of  scatterers,  i.e/Xthe  quasi¬ 
entire  function,  may  find  concise  description  in  a  fonsVimilar  to 
extended  physical  optics,  wherein  the  shadowing  effects  inherent  in  the 
standard  definition  of  physical  optics  are  ignored. 
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IOTWDDUCTION 


Over  the  years  since  the  invention  of  the  singularity  expansion 
method  (SEN),1  there  has  been  some  uncertainty  surrounding  the  identity 
and  role  of  the  "possible  entire  function"  term  in  the  singularity  ex¬ 
pansion  of  a  scattered s  response.  The  fact  that  the  entire  function 
does  not  have  a  unique  definition  causes  the  role  of  this  component  to 
be  somewhat  vague.  The  results  presented  herein  represent  an  attempt  to 
grasp  the  meaning  and  utility  of  the  entire  function.  The  present  in¬ 
vestigation  seeks  to  answer  the  following  questions:  (1)  Is  fun  entire 
function  component  required  for  the  accurate  representation  of  transient 
scattering  responses?  (2)  What  is  the  physical  significance  of  an 
entire  function  component?  (3)  For  applications  such  as  resonance  esti¬ 
mation  and  radar  target  identification,  how  can  an  entire  function  com¬ 
ponent  be  represented,  dealt  with,  or  effectively  used? 

Entire  functions  are  frequently  found  in  the  frequency-domain 
responses  of  circuits  where  they  represent  forced  components  in  the 
responses.  Forced  components  are  distinguished  from  the  unforced  or 
natural  responses  of  a  system  in  that  they  exist  only  while  the  system 
is  under  excitation,  whereas  the  unforced  response  may  exist  independ¬ 
ently  of  the  excitation.  Although  the  unforced  response  can  exist  inde¬ 
pendently,  in  all  practical  situations  it  is  the  result  of  some  past  (or 
present)  excitation.  Therefore,  the  label  "unforced"  is  misleading.  In 
fact,  the  forced  response  can  be  thought  of  as  that  part  of  the  total 
response  which  quickly  damps  out  after  the  excitation  is  removed  and  the 
unforced  response  is  the  remainder  with  little  real  distinction  between 
the  two  conponents .  The  recognition  of  this  similarity  between  the 
forced  and  unforced  components  is  crucial  for  understanding  one  facet  of 
the  entire  function's  role  in  transient  scattering. 

As  is  shown  in  Sec.  2  for  the  far-field  response  of  a  conducting 
sphere  under  plane  wave  illumination,  an  entire  function  component  is 
not  required.  On  the  other  hand,  if  the  sphere  is  deformed  through  the 
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prolate  spheroids  into  an  infinitely  thin  filament,  an  entire  function 
is  required.  As  the  poles  of  the  sphere  are  tracked  from  the  sphere  to 
the  filament,  some  of  the  poles  approach  the  imaginary  axis  in  the 
s-plane  while  the  rest  tend  to  s  -  -®  .  The  poles  that  withdraw  to 
infinity  give  rise  to  an  entire  function  whicn  can  be  thought  of  as  the 
forced  component  of  the  filament's  response.  In  the  time-domain,  this 
forced  component  is  a  square-pulse  function,  i.e.,  a  time-limited  func¬ 
tion,  that  is  non-zero  only  while  the  incident  plane  wave  is  on  the 
scatterer. 

In  Sec.  3,  the  approximation  of  time-limited  functions  with  sums 
of  highly  damped  exponentials  is  found  to  be  theoretically  feasible 
although  in  practice  such  approximation  encounters  numerical  difficul¬ 
ties.  In  the  frequency-domain,  this  type  of  approximation  corresponds 
to  representing  an  entire  function  by  a  pole  series  involving  poles  with 
large  negative  real  parts.  Hie  practical  difficulties  inherent  with 
such  a  representation  points  out  the  effectiveness  of  the  entire  func¬ 
tion  as  an  approximation.  For  distributed  systems,  the  entire  function 
component  can  be  used  as  an  approximate  representation  for  the  contribu¬ 
tion  of  an  infinitely  large  group  of  poles  that  recede  to  infinity  in 
the  complex  frequency  plane.  Hie  forced  components  in  the  responses  of 
circuits  are,  in  fact,  based  on  such  approximations.  In  an  absolute 
sense,  an  entire  function  may  never  be  required  in  the  responses  of  real 
systems.  In  a  practical  sense,  however,  the  entire  function  can  provide 
a  welcome  refuge  from  all  the  numerical  problems  posed  by  attempting  to 
represent  transient  responses  by  pole  series  with  infinite  numbers  of 
terms. 

Such  practical  considerations  motivate  the  present  investigation. 
Radar  target  identification  has  been  suggested  as  a  possible  application 
of  the  SEN  representation  wherein  the  poles  of  a  radar  scatterer  are 
used  to  identify  it.  Poles  are  particularly  useful  features  for  target 
identification  because  they  are  invariant  with  aspect  angle.  Hie 


authors  first  became  acquainted  with  the  question  of  the  entire  function 
while  estimating  resonances  from  measurements  of  transient  scattering 
from  thin  cylinders.  We  noticed  that  only  when  the  early-time  portions 
of  the  back  scattered  returns  were  excluded  from  the  resonance  estima¬ 
tion  procedure  (Prony's  method)  could  consistent  estimates  of  the  reson¬ 
ances  be  achieved.  At  the  time,  we  speculated  that  the  early-time  data 
was  not  composed  purely  of  the  unforced  response,  but  also  included  a 
forced  ccmponent  which  could  not  be  adequately  modeled  by  the  finite 
pole  series  being  matched  to  the  data.  The  results  of  the  present  in¬ 
vestigation  confirm  this  initial  speculation. 

Because  an  infinite  number  of  degrees  of  freedom  exist  in  the 
residues  of  a  scatterer's  pole  series,  there  is  considerable  freedom  in 
the  choice  of  the  entire  function.  The  residues  can,  in  many  cases,  be 
adjusted  to  absorb  or  produce  an  additional  arbitrary  entire  function  in 
the  representation.  In  Sec.  3,  it  is  shown  that  the  Laplace  transform 
of  a  time-limited  function  can  be  represented  to  any  degree  of  accuracy 
by  a  pole  series  formed  from  poles  with  large  negative  real  parts. 

Since  most  scatterers  possess  an  infinite  number  of  such  poles,  any 
entire  function  resulting  from  an  arbitrary  time-limited  function  in  the 
early-time  response  can  be  expanded  in  those  poles  to  any  degree  of 
accuracy.  As  a  result,  there  is  usually  a  variety  of  choices  for  entire 
function  and  one  simply  has  to  make  a  good  choice. 

2 

Heyman  and  Felsen  have  suggested  that  for  many  types  of  scat¬ 
terers,  such  as  the  sphere,  pole  series  are  not  very  useful  in  the 
representation  of  the  early-time  response.  Instead,  they  suggest  that, 
for  many  scatterers,  the  responses  can  be  more  efficiently  represented 
if  groups  of  resonances  are  combined  into  the  wavefronts  and  rays  of  the 
geometrical  theory  of  diffraction  for  the  early-time  responses,  while 
the  SEM  pole  series  can  be  retained  for  late  times  if  the  late-time 


response  is  significant.  Wavefronts  are  nothing  more  than  entire  func¬ 
tion  components .  Specifically  for  the  conducting  sphere,  Uberall  and 
Gaunaurd^  have  shown  hew  the  scattered  fields  of  the  conducting  sphere 
can  be  split  into  a  geometrical -optics  and  creeping-wave  components. 

The  geometrical-optics  component  corresponds  to  an  entire  function  and 
the  creeping  waves  to  a  pole  series  in  the  frequency  domain.  Mathemat¬ 
ically,  this  entire  function  has  no  meaning  since  it  could  be  excluded 
from  the  representation  altogether.  Numerically,  the  inclusion  of  this 
component  accelerates  the  convergence  of  the  representation.  Physic¬ 
ally,  this  component  can  be  given  meaningful  interpretation  as  the  re¬ 
flection  from  the  apex  of  the  sphere.  This  choice  for  the  entire  func¬ 
tion  is,  of  course,  one  among  many. 

4  5 

Morgan  and  Pearson  have  suggested  that  the  entire  function  is 

required  in  addition  to  the  pole  series  to  accurately  describe  the  scat- 

terer's  response.  The  results  of  the  present  investigation  do  not  bear 

out  the  absolute  necessity  of  the  entire  function  in  the  responses  of 

scatterers  but  only  its  usefulness  in  the  accurate  and  efficient  repre- 

4 

sentation  of  the  early-time  response.  Morgan's  entire  function  repre¬ 
sents  another  possibly  good  choice  for  the  entire  function. 

For  many  scatterers,  it  is  possible  to  distinguish  two  groups  of 
poles:  the  highly-damped  poles  and  the  dominant  poles.  The  contribu¬ 
tions  of  the  highly-damped  poles  can  be  grouped  into  what  we  will  call  a 
"quasi-entire  function"  component.  In  Sec.  4,  the  quasi-entire  compo¬ 
nents  of  two  scatterers,  a  conducting  cylinder  and  loop,  are  numerically 
derived  by  subtracting  the  dominant  pole  series  from  the  total  re¬ 
sponses.  Both  class  1  and  class  2  coupling  coefficients  are  used  to 
form  the  pole  series  which  results  in  two  definitions  of  the  quasi¬ 
entire  function  for  each  case,  thus  illustrating  the  non-uniqueness  of 
this  response  component.  The  quasi-entire  function  may  be  regarded  as  a 
forced  response  component  in  the  same  spirit  as  the  forced  responses  of 

circuits.  As  the  results  of  Sec.  4  illustrate,  this  forced  component 

4  6 

resembles  the  extended  physical  optics  '  response  of  the  scatterer. 
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EXAMPLES  AND  COUNTEREXAMPLES  OF  REQUIRED  ENTIRE  FUNCTION  COMPO¬ 
NENTS  IN  TRANSIENT  RESPONSES 

To  illustrate  the  basic  issues  under  consideration,  some  simple 
systems  that  can  be  analyzed  by  exact  methods  allow  some  preliminary 
conclusions  to  be  made.  First,  a  circuit  shown  in  Fig.  2.1,  whose 
transfer  function  contains  both  an  entire  function  component  and  a  reso¬ 
nant  component  is  considered.  The  transfer  function  of  this  circuit  can 
be  written  as 


Vin(8)  -1  -1 

H(s)  -  y— '-fsy  -  1  ♦  Ajfs  -  sx)  +  AjU  -  s2) 


(2.1) 


where 


l,2  "  J 
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*  L  S2  “ 


R 


Figure  2.1.  A  network  whose  transfer  function  contains  an  entire 
function  and  a  pole  series. 
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with  1/(LC)  >  fV(2L)^  .  The  first  term  of  Eq.  2.1  is  a  (constant) 
entire  function  of  the  Laplace  transform  variable  s  -  <r  +  j«  .  The 
presence  of  this  term  causes  the  output  of  the  circuit  to  contain  the 
familiar  forced  component .  The  remaining  two  terms  comprise  a  finite 
pole  series  that  constitutes  the  unforced  or  natural  response  of  the 
circuit.  The  entire  function  in  this  case  is  required  for  the  accurate 
representation  of  the  response. 

Next,  consider  the  direct  connection  between  input  and  output 
formed  by  replacing  the  R  element  of  Fig.  2.1  by  a  short  and  the  L 
and  C  elements  with  an  open.  The  transfer  function  of  this  network  is 
H(s)  -  1  and  the  response  consists  only  of  a  forced  entire-function 
component.  This  direct  connection  can  be  thought  of  as  transmission 
line  observed  at  low  frequencies.  At  high  frequencies,  a  transmission 
line  of  length  L  with  a  voltage  generator,  vin  ,  at  the  input  and  an 
open  at  the  output  as  shown  in  Fig.  2.2  has  a  transfer  function  given  by 

H(s)  "  cos ( — jsL/c )  (2,2) 

where  c  is  the  speed  of  light.  The  transfer  function  has  poles, 


Figure  2.2.  A  transmission  line  as  a  network  whose  transfer  function 
can  be  represented  exactly  as  an  infinite  pole  series 
alone  or  approximately  at  low  frequencies  by  an  entire 
function  alone. 
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,  j  (2n  +  l)cn 
sn  "  3  Zl 

and  can  be  expanded  into  a  pure  pole  series  with  an  infinite  number  of 
terms.  Note  that  at  low  frequencies  Eq.  2.2  can  be  approximated  as 
H(s)  ■  1  .  It  then  becomes  obvious  that  the  entire  function  found  in 
circuit  theory  is  a  result  of  the  low  frequency  approximations  that  form 
the  basis  of  circuit  theory.  A  large  number  of  high  frequency  resonant 
contributions  can  be  combined  and  approximately  represented  as  an  entire 
function  accurately  within  limited  region  of  the  complex-frequency 
plane. 


The  above  considerations  indicate  that  an  entire  function  contri¬ 
bution  may  never  exist  in  reality  but  is  merely  a  convenient  approximate 
representation  for  a  large  number  of  high-frequency  resonant  contribu¬ 
tions.  The  question  to  answer  then  becomes:  how  is  the  entire  function 
to  be  defined  in  specific  cases?  The  answer  to  this  question  appears  to 
be  that  no  precise  or  unique  definition  is  ever  possible.  However ,  in 
certain  extreme  system  configurations,  the  system  poles  will  frequently 
cluster  into  widely  separated  regions  of  the  complex  plane,  in  such 
cases,  the  definition  of  the  entire  function  can  take  on  a  measure  of 
precision  and  uniqueness.  One  such  extreme  configuration  is  the  infi¬ 
nitely  thin  wire  scatterer  which  is  explored  below.  Other  such  config¬ 
urations  are  circuits  of  which  an  example  has  already  been  given.  Next, 
we  will  consider  two  scatterer  geometries,  namely  the  sphere  and  the 
infinitely  thin  wire,  which  illustrate  the  extremes  with  regard  to  the 
entire  function  question. 

The  conducting  sphere  is  one  scatterer  that  can  be  solved  entirely 
by  analytical  techniques.  Moreover,  the  poles  of  the  sphere  are  found 
by  finding  the  zeros  of  polynomials,  as  is  done  for  circuits,  instead  of 
more  complicated  functions  as  is  the  case  for  most,  if  not  all  other 
scatterers.  The  scattered  far  field  of  a  perfectly  conducting  sphere  in 
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free  space  with  an  +x-polarized  incident  plane  wave  traveling  in  the  +z 

1  7— Q 

direction  can  be  written  as' 


Is  -  _  j  exp(-jkR)  2n  +  1 
J  R  W  n(n  +  1 


n-1 


T7 


’ 

jn(ka) 

P*(COS  0) 

kh^2,(ka) 

sin  6 

co,  4  -  4 


V 

lkajn(ka)]' 

dP*(cos  0) 

k[kah£2,(ka)]' 

d0 

COS 
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pJ(COS 

sin  6 


0) 


sin 


(2.3) 


where  p*  is  the  associated  Legendre  function  of  the  first  kind  of 


order  1  and  degree  n  ,  jn  and  h^  '  are  spherical  Bessel  functions 
of  the  first  and  third  kinds,  respectively,  a  is  the  radius  of  the 
sphere,  k  «  <n/c  -  -js/c  is  the  wave  number,  and  R  is  the  distance  to 
the  far  field  point.  The  derivatives  of  the  spherical  Bessel  functions 
appearing  in  Eg.  2.3  are  defined  by 


lkazn(ka)]'  -  ^  lxzn(x)]  x_k#  -  ^(ka)  -  ^  zn(ka) 


where  zR  denotes  the  Bessel  function.  Equation  2.3  can  be  converted 
into  a  pole  series  plus  an  entire  function  by  expanding  the  factors 


lkajn(ka)]' 
k[kah£2,(ka)]' 

into  pole  series  and  entire  function  terms  where  the  spherical  Bessel 
functions  are  defined  by 

h^’u)  -  j-"-1  £SEm  £  ai(-2iz)-i 

i-0 

h<2l(Z)  -  j-1  asking  ai(2j2ri 

i-0 


1  + 


[kahJ^Oca)]' 

[kah£2,(ka)]' 


(2.4b) 


and 


ai 


(n  +  i) l 


ii(n  -  i)l 


For  each  n  ,  Eqs.  2.4a  and  2.4b  can  be  expanded  into  a  pole  series  plus 
the  entire  function,  [1  -  exp( j2ka) ]/(2k)  .  These  entire  functions  can 
be  sunned  over  n  to  give  a  total  entire  function  component  in  the 
sphere's  response.  Upon  examination,  this  component  displays  some 
rather  peculiar  features,  has  no  apparent  physical  interpretation,  and 
is  regarded  as  a  mathematical  artifact. 

Although  when  each  of  the  factors  in  Eq.  2.4  are  taken  individu¬ 
ally,  an  entire  function  term  must  exist  for  each  term  of  Eq.  2.3;  this 
does  not  indicate  that  the  total  response  requires  an  entire  function. 

In  Sec.  3,  it  is  shown  that  the  residues  of  a  pole  series  with  a  large 
number  of  deeply-imbedded  poles  can  always  be  adjusted  to  absorb  the 
entire  function  into  the  pole  series.  Since  the  sphere  possesses  such  a 
set  of  poles,  the  response  of  the  sphere  should  be  representable  to  any 


2-5 


degree  of  accuracy  by  a  pole  series  alone.  It  is  always  possible, 
though,  to  rearrange  the  equations  so  that  an  entire  function  component 
is  retained.  In  this  case,  the  residues  of  the  resulting  pole  series 
are  altered  to  account  for  the  additional  component . 

Next  the  infinitely  thin  conducting  wire  or  filament  scatterer  is 
considered.  This  scatterer  requires  an  entire  function  component,  in 
addition  to  a  constant  coefficient  pole  series,  to  accurately  represent 
its  scattered  field  and  thus  it  provides  a  counter  example  to  that  of 
the  sphere.  Finally,  the  transition  from  sphere  through  the  prolate 
spheroids  to  the  filament  will  be  explored  and  there  the  specific  causes 
of  the  difference  between  the  sphere  and  the  filament  with  regard  to  the 
entire  function  will  become  clear. 

To  develop  a  solution  for  the  axial  currents  on  a  filament  with  a 
plane  wave  incident  field,  we  start  with  an  asymptotic  solution  obtained 
for  an  infinitely  long  conducting  cylinder  driven  by  voltage  across  an 
infinitely  thin  gap,  given  by  Chen  and  Keller.10  This  solution,  which 
is  valid  for  large  k|z  -  z'|  ,  can  be  written  as  an  admittance  which, 
when  multiplied  by  the  gap  voltage,  gives  the  current: 

.is - .«P<-8|2  -  2- I/O  -  (2.5) 

^  log(2j jz  -  z'|/a)  -  log(-jI^ sa/c) 

where  h  -  is  the  impedance  of  free  space  in  which  the  cylinder 

is  immersed,  a  is  the  radius  of  the  cylinder,  r  »  1.781  is  the  expo¬ 
nential  of  Euler's  constant,  z'  is  the  z  coordinate  of  the  gap  and 
z  is  the  coordinate  for  the  location  of  the  current.  The  tilde  denotes 
a  Laplace  transformed  quantity. 

For  very  small  ka  ,  Eq.  2.5  can  be  written  in  a  simpler  form: 

mm 

Y-(z,z',s)  «  Yq  exp(-s|z  -  z'|/c)  (2.6) 
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where 


Y  -  -  2n 
o  “  n  log  wQ a/c  ' 

and  «0  is  a  constant  such  that  wQ  a  w  .  in  transient  analysis,  w 
will  vary  over  sene  range.  In  this  case,  «Q  is  set  to  the  middle  of 
this  range. 

Equation  2.6  is  valid  even  for  k|z  -  z' |  «  1  when  ka  is  very 
small.  However,  in  a  region  very  close  to  the  gap,  Eq.  2.6  does  not 
adequately  account  for  the  logarithmic  singularity  in  the  current.  Near 
the  gap,  the  current  is  given  by1® 

5“  «  -  log  k|z  -  z'  | 

As  ka  -»  0  ,  this  component  has  a  vanishingly  small  effect  compared  to 
the  other  components  of  current  given  by  Eq.  2.6  and,  hence,  it  can  be 
neglected. 

The  solutions  of  Eqs.  2.5  and  2.6  assume  that  a  gap,  i.e.,  an 
opening  in  the  metal,  exists,  and  that  the  gap  voltage  has  no  variation 
in  4  so  that  no  ^-directed  currents  are  excited  on  the  cylinder.  Some 
resolvable  differences  then  exist  between  the  gapless  filament  under 
excitation  by  an  arbitrary  plane  wave  and  the  cylinder.  The  gap  in  the 
cylinder  can  be  closed  without  any  change  in  the  exterior  fields.  This 
is  done  by  defining  fields  of  the  gap-excited  cylinder  as  scattered 
fields,  closing  the  gap  with  a  perfect  conductor,  and  introducing  an 
incident  field  to  cancel  with  the  electric  field  corresponding  to  the 
gap  voltage  on  the  perfect  conductor.  It  will  be  assumed  that  the  4 
components  of  the  incident  plane  wave  electric  field  on  the  cylinder's 
surface  and  the  ^-directed  currents  produced  by  these  components  have 
negligible  effects  on  the  overall  response  of  the  filament.  This  is  a 
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standard  thin-wire  assumption  which  is  known  to  produce  correct  results 
as  ka  approaches  zero.  The  ♦-variation  of  the  incident  field  obvi¬ 
ously  becomes  negligible  as  ka  approaches  zero. 

Applying  the  inverse  Laplace  transform  to  Eg.  2.6  gives  the  time 
domain  form: 


Y*(z,z',t)  «  Y0S(t  -  |z  -  z'|/c)  (2.7) 

where  i(t)  is  the  Dirac  delta  function.  Equation  2.7  approaches  the 
exact  solution  away  from  the  gap  for  the  time-domain  current  on  an 
infinitely  long  filament  excited  locally  at  z'  by  a  delta  function 
incidence  field.  Equation  2.7  represents  two  current  impulses  traveling 
in  both  directions  away  from  the  excitation  point. 

Equation  2.7  can  be  used  to  determine  the  currents  on  a  filament 
of  finite  length.  The  reflections  of  the  current  pulses  off  the  ends  of 
a  finite  filament  can  be  modeled  with  additional  sources  (infinitely 
thin  gaps  with  voltages  across  them)  situated  at  the  ends  and  defined 
such  that  the  currents  at  the  ends  are  always  zero.  Chen  and  Keller^ 
give  a  detailed  description  of  the  modeling  of  the  ends  of  the  cylinder. 
The  time-domain  current  on  a  filament  of  length  L  excited  by  a  local¬ 
ized  E-field  or  voltage  source  at  z'  can  then  be  written  as 

• 

Y(z,z',t)  -  Yq  ^  js(ct  -  z  +  z'-  2nL)  +  8(ct  +  z  -  z*  -  2nL) 
n-0 

-  4{ct  -  z  -  z'  -  (2n  +  1)L] 

-  8[ct  +  z  +  z'  -  (2n  +  l)L)^u(t)  (2.8) 


2-8 


where  u(t)  is  the  unit  step  function.  Equation  2.8  represents  the  two 
current  inpulses  traveling  away  from  z'  and  recurrently  reflecting  off 
the  ends  of  the  finite  filament. 


Filament  currents  due  to  an  incident  impulsive  plane  wave,  as 
shown  in  Fig.  2.3,  can  be  found  by  superposition,  using 


i  Y2 

I(z,t)  -  E*  /  Y(z,z’ ,t  +  z'  cos  0')  dz'  (2.9) 

-L/2 


where  is  the  z  cooponent  of  the  incident  electric  field.  Simi¬ 
larly,  the  scattered  far  field  can  be  found  from 

L/2 

E0(e,t)  -  ^  J  I{z,t  +  z  cos  0  -  R)  dz  (2.10) 

-L/2 

where  0  is  the  scattering  angle  and  R  is  the  distance  to  the  far 
field.  Hie  currents  and  fields  given  by  Eqs.  2.9  and  2.10  are  composed 
of  the  superposition  of  a  number  of  step  functions  in  time. 


L 


* 


Figure  2.3.  Geometry  for  scattering  by  filament  scatterer  of  length  L. 


Figure  2.4  illustrates  the  time  history  of  the  backscattered  far 
field  produced  by  an  infinitely  thin  filament  of  length  L  with  impul¬ 
sive  plane  wave  excitation  as  found  using  Eq.  2.10.  In  Fig.  2.4,  the 
responses  of  a  perfectly  conducting  sphere-capped  cylinder  with  a 
length-to-diameter  ratio  of  30  are  provided  for  comparison  to  those  of 
the  filament.  The  cylinders  response's  were  obtained  using  a  moment- 
method  solution  of  the  electric  field  integral  equation.  Note  that  all 
of  the  features,  except  for  the  damping  and  amplitudes,  of  the  fila¬ 
ment's  and  cylinder's  responses  agree  well.  The  late-time  portions  of 
the  filament  responses  are  completely  periodic,  whereas  in  the  early- 
time  an  aperiodicity  is  noted.  The  periodic  late-time  portion  of  the 
response  can  be  expanded  in  a  Fourier  series  which  is  nothing  more  than 
the  pole  series  employing  the  dominant  filament  resonances.  However, 
this  series  cannot  model  the  aperiodicities  in  the  early-time  portions 
of  the  filament  responses.  If  the  periodic  late-time  response  is  con¬ 
tinued  into  the  early-time  and  subtracted  from  the  total  response  in  the 
early-time,  the  residual  is  a  time-limited  function  whose  Laplace  trans¬ 
form  is  an  entire  function. 

In  contrast  to  the  sphere,  a  filament  scatterer  requires  an  entire 
function  component  in  its  far  field  response  to  plane  wave  excitation. 
What  is  the  difference  between  the  two  cases?  As  is  shown  below,  some 
of  the  sphere's  poles  withdraw  to  s  -  -•  as  the  sphere  is  deformed 
through  prolate  spheroids  toward  the  infinitely  thin  filament.  It  is 
these  large  frequency  poles  whose  contributions  may  be  approximated 
quite  closely  by  an  entire  function  in  the  manner  discussed  above.  In 
the  limit  as  ka  -»  0  ,  the  approximation  becomes  exact. 

Figures  2.5  and  2.6  show  selected  pole  loci  for  the  prolate  spher¬ 
oid  as  it  varies  from  the  sphere  to  the  filament  scatterer.  The  geom¬ 
etry  of  the  prolate  spheroid  is  indicated  in  Fig.  2.7.  In  cylindrical 

2  2  2  2 

coordinates,  the  equation  for  the  spheroid  is  p  /a  +  z/b  -  1  .  The 
results  for  b/a  <  •  were  obtained  with  BRFD11  which  is  a  scattering 
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(b)  Sphere-capped  cylinder  with  a  length-to-diameter  ratio  of  30 
(L/D  -  30),  9  -  90# 

Figure  2.4.  Time  history  of  backscattered  far  fields  due  to  impulsive 
plane  wave  excitation. 


Figure  2.5.  Prolate  spheroid  pole  loci  for  surface  currents  nodes 
which  are  constant  in  +. 
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Figure  2.6.  Prolate  spheroid  pole  loci  for  surface  currents  modes 
which  have  first  harmonic  variation  in  ♦,  i.e.,  the 
modes  vary  as  exp(+j$). 


code  for  general  conducting  bodies  of  revolution  and  is  based  on  a 
method-of-moments  solution  of  the  frequency-domain  electric  field  inte¬ 
gral  equation.  These  results  indicate  that  certain  poles  of  the  sphere, 
namely  those  poles  in  the  layer  nearest  the  imaginary  axis  which  are 
associated  with  current  modes  that  are  constant  in  4  ,  migrate  to  the 
imaginary  axis  in  the  s-plane  as  the  sphere  is  deformed  continuously 
through  the  prolate  spheroids  to  arrive  finally  at  the  filament  scat¬ 
tered  All  other  poles  withdraw  to  s  -  .  The  periodic  components 

of  the  filament  responses  shown  in  Fig.  2.4  must  be  composed  of  pole 
series  of  the  filament  which  is  nothing  more  than  a  Fourier  series.  The 
early-time  aperiodic  component  must  be  due  to  those  poles  that  have 
withdrawn  to  s  -  -•  .  Hence,  the  poles  that  flee  to  infinity  must  give 
rise  to  the  entire  function  in  the  filament's  frequency-domain  response. 

When  the  prolate  spheroid  is  thin,  but  not  infinitely  thin,  there 
exists  a  group  of  poles  which  are  deeply- imbedded  in  the  left  half  of 
the  s-plane.  These  poles  give  rise  to  a  quasi-entire  function  in  the 
scattered s  frequency-domain  response  and  are  distinguished  from  the 
dominant  poles  which  lie  near  the  imaginary  axis.  For  general  scatterer 
geometries,  it  should  always  be  possible  to  define  a  quasi-entire  func¬ 
tion  whenever  a  set  of  dominant  poles  are  clearly  distinguishable  from 
the  rest.  The  numerically-derived  quasi-entire  functions  of  selected 
thin  scatterers  are  examined  in  Sec.  4. 
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3  APPROXIMATE  REPRESENTATION  OF  TIME-LIMITED  FUNCTIONS  BY  EXPONEN¬ 
TIAL  SUMS  ~  ~ 

In  Sec.  2,  we  have  seen  that  the  contribution  of  a  group  of  large 
frequency  poles  is  an  effectively  time-limited  component  in  a  scat¬ 
tered  s  response  which  corresponds  to  a  quasi-entire  function.  In  this 
section,  the  representation  of  actual  time-limited  functions  by  exponen¬ 
tial  series,  which  are  the  time-domain  representatives  of  pole  series, 
is  considered.  It  is  found  that  time-limited  functions  can  be  repre¬ 
sented  to  any  degree  of  accuracy  with  a  exponential  series  by  including 
a  sufficient  number  of  terms  and  that  the  best  approximations  are 
obtained  if  poles  that  are  deeply-imbedded  in  the  left  half  of  the 
s-plane  (Re  s  large  and  negative)  are  employed  in  the  exponential 
series.  Hie  fact  that  the  pole  sets  of  scatterers  usually  include  many 
such  deeply-imbedded  poles  suggests  that  the  residues  of  the  complete 
pole  series  for  the  scattered s  response  can  always  be  adjusted  to 
incorporate  any  time-limited,  forced  component  that  may  exist. 

First  consider  the  approximation  of  the  exponential  entire  func¬ 
tion  by  a  pole  series.  As  is  shown  in  Appendix  A, 

N 

exp(z)  a  if1  7T  <zi  ~  z>~1  (3.1) 

i-1 


where  the  following  conditions  hold 


|  Z  |  «  N 


| zi  -  N|  «  N  for  all  i 


1 
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t 

i-1 


(z*  -  N) 


«  1 


(3.2a) 

(3.2b) 


(3.2c) 
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Note  that  Eg.  3.1  can  be  expanded  in  partial  fractions  to  form  a  pole 
series.  The  z.  can  have  any  distribution,  so  long  as  the  conditions 
in  Egs.  3.2b  and  3.2c  are  satisfied,  lhus,  in  the  finite  complex  plane, 
the  exponential  entire  function  can  be  approximated  to  any  degree  of 
accuracy  by  a  pole  series  with  a  sufficient  number  of  terms.  Conditions 
in  Eqs.  3.2b  and  3.2c  are  roughly  equivalent  to  the  constraint  that 
z^  =  N. 

The  Laplace  transform, 
b 

F(s)  -  /  f (t)e“st  dt  ,  0  <  a  <  b 

a 

of  the  time-limited  function,  f(t}  ,  can  be  approximated  to  any  degree 
of  accuracy  by  choosing  M  sufficiently  large  in  the  sum 

M 

F(s)  a  f(tk)  exp(-stk)At  (3.3) 

k-1 

where 

tk  -  •  +(*  -  j)“ 


and 


At  -  (b  -  a)/M 

Here  we  assume  that  f(t)  satisfies  the  usual  conditions  for  the  exis¬ 
tence  of  the  transform.  Replacing  z  with  -stk  and  z^  with 
z^k  -  — siktk  in  Eq.  3.1,  and  substituting  into  Eq.  3.3  gives 
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M  N 

F(s)  a  f(tk)AtNN7T  <zik  "  Z)~X  (3.4) 

k-1  i-1 


where  the  zik  are  chosen  so  that  the  conditions  in  Eqs.  3.2a  to  3.2c 
are  satisfied.  Equation  3.4  can  be  expanded  into  partial  fractions  to 
provide  a  pole  series  for  F(s)  .  Therefore,  the  following  has  been 
shown:  The  Laplace  transform  of  any  time-limited  function  can  be 
approximated  to  any  degree  of  accuracy  over  a  finite  region  of  the  com¬ 
plex  plane  with  a  pole  series.  The  restriction  to  the  finite  region  of 
the  complex  plane  is  equivalent  to  the  restriction  to  band-limited 
cases. 


If  Eq.  3.4  is  to  be  a  good  approximation,  the  s^  should  satisfy 
sik  “ 

for  all  i  and  k  which  is  equivalent  to  saying  that  the  poles  should 
be  deeply  imbedded  in  the  left  half  of  the  s-plane.  In  general,  scat- 
terers  possess  an  infinite  number  of  poles  that  satisfy  this 
requirement. 

The  above  observations  indicate  that,  at  least  for  bandlimited 
cases,  the  residues  or  coefficients  of  the  pole  series  for  a  scattered s 
response  to  an  impulsive  plane  wave  can  always  be  adjusted  to  incorpo¬ 
rate  any  time-limited,  forced  component  that  may  exist  while  the  inci¬ 
dent  plane  wave  is  on  the  scatterer,  and  thus  that  the  natural  modes 
form  a  "complete"  basis  for  the  scatterer 's  responses,  unfortunately, 
this  completeness  property  of  the  natural  modes  cannot  be  demonstrated 
with  absolute  rigor  since  the  above  argument  is  based  on  an  assumed 
rough  knowledge  of  the  pole  distribution  in  the  complex  frequency  plane. 
For  general  scatterers,  the  authors'  knowledge  is  quite  limited  in  this 
regard. 
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To  more  concretely  explore  the  question  of  the  representation  of 
time-limited  entire  functions  as  discussed  above,  the  following  numeri¬ 
cal  experiment  was  performed.  A  pole  series  of  the  form 

N 

F(s)  Ai(i  -  l)l(s  +  N)1  (3.5) 

i-1 

which  in  the  time  domain  is 

N 

f(t)  -  2  Aitl_1  exP(-Nt>  (3.6) 

i-1 

is  fit  to  the  pulse  function  u(t)  -  u(t  -  1)  in  a  least-squares  sense 
over  the  interval  (0,-)  by  adjusting  the  Ai  ,  as  described  in  Appendix 
B.  In  Eq.  3.5,  F(s)  has  an  nth  order  deeply-imbedded  pole  at 
s  «  -N  .  Figure  3.1  shows  the  results  for  N  -  2  ,  4,  and  6.  This  type 
of  approximation  on  an  infinite  interval  cannot  be  achieved  with  a 
Fourier  series,  which  can  only  represent  periodic  functions,  except  in 
the  limit  as  the  Fourier  series  becomes  the  Fourier  transform.  The 
Fourier  series  is  equivalent  to  a  pole  series  whose  poles  are  evenly 
spaced  along  the  imaginary  axis  in  the  s-plane.  Figure  3.1  illustrates 
the  use  of  deeply-imbedded  poles  in  a  pole  series  for  the  approximation 
of  an  aperiodic,  time-limited  function  on  the  infinite  interval.  The 
practicality  of  using  such  a  representation,  though,  is  questionable 
since  as  the  poles  move  deeper  into  the  left  half  of  the  s-plane,  one 
encounters  numerical  problems.  One  way  to  avoid  these  problems  is  to 
abandon  the  pole  series  and  resort  to  the  entire  function  representation 
obtained  by  taking  the  Laplace  transform  of  the  pulse  function.  There¬ 
fore,  as  a  practical  matter,  pole  series  using  deeply-imbedded  poles  are 
abandoned  in  favor  of  other  representations,  one  of  which  is  the  entire- 
function  representation. 
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Figure  3 


> 


(a)  Sum  of  2  exponentials 


1.  Approximation  of  pulse  function  by  a  sum  of  highly-damped 
exponentials.  The  dashed  lines  show  the  pulse  function. 
The  solid  lines  show  the  approximation. 
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(c)  Sum  of  6  exponentials 
Figure  3.1  (Cont.) 


Ordinarily,  one  may  resort  to  other  representations  for  the  quasi¬ 
entire  function  with  impunity  except  while  attempting  to  estimate  the 

12 

resonances  with  techniques  like  Prony's  method.  For  resonance  estima¬ 
tion,  it  may  be  feasible  to  assign  a  number  of  highly-damped  poles  to 
account  for  the  quasi-entire  function  in  a  constrained  resonance  estima¬ 
tion  procedure  like  the  constrained  version  of  Prony's  method  described 
in  Appendix  C. 
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NUMERICAL  EXAMPLES  OF  TOE  QUASI -ENTIRE  FUNCTION 
In  this  section,  the  numerically-derived  quasi-entire  function, 
which  was  discussed  in  Sec.  1,  of  the  conducting  loop  (torus)  and  the 
sphere-capped  cylinder  are  examined.  Both  of  these  scatterers  possess  a 
clearly  distinguishable  set  of  dominant  resonances  when  they  are  suffi¬ 
ciently  thin. 

The  method  for  isolating  the  quasi-entire  function  consists  of 
subtracting  the  pole  series  constructed  with  the  dominant  poles  from  the 
total  response.  The  pole  series  is  constructed  through  use  of  the  elec¬ 
tric  field  integral  equation  (EFIE)  for  a  perfectly-conducting  body, 
which  can  be  written  as 

_  A  . 

rj  -  n  X  E1  (4.1) 

where  J  is  the  surface  current  density,  I1  is  the  incident  electric 
A 

field,  n  is  the  outward-pointing  unit  vector  normal  to  the  scatterer's 
surface,  and  r  is  a  linear  operator  defined  by 

13  -  A  x 

The  magnetic  vector  potential,  A  ,  and  electric  scalar  potential,  #  , 
are  defined  as 

A(r)  -  f  J(r')  e3q?4^kR)  dS' 

S 

*(r)  -  J V'-J(r')  dS' 

S 
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where  the  integrals  are  over  the  scattered s  surface,  S  ,  and 
R-  |r  -  r'  | . 


The  method  of  moments1^  is  allied  to  convert  the  EFIE  into  a  set 
of  matrix  equations.  The  conversion  is  accomplished  by  first  expanding 
the  current  on  a  set  of  basis  functions,  ^ffl(r),  m  -  1,  . ,  as 


J(I'>  ~Yj  (4-2) 

m 

Next,  define  an  inner  product 

(f,g)  -  f I-g  dS 
S 

Finally,  a  set  of  testing  functions  ^g^  m  -  1,  . is  defined  and 
a  set  of  equations  is  formed  by  taking  the  inner  product  of  Eq.  4.1  with 
each  of  the  testing  functions.  Using  Eq.  4.2,  these  equations  can  be 
written  as 


M 

12  Jm(rIm'^n)  "  n  "  1'2,...,M 

unI 

These  equations  can  be  written  in  matrix  form  as 
ZI  -  V 

where  the  nth  element  of  V  is  (A  x  E1,  gfl),  the  mth  element  of  I  is 

J  ,  and  the  nth-mth  element  of  Z  is  (rf„,g  ). 
m  —  —  m  n 
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Assuming  only  first  order  dominant  poles,  the  response  for  the 

14  15 

scatterer  currents  may  be  written  in  the  form  ' 


I(s)  -  (s  -  s  )  \  h  (s)  +  W(s) 

bam  01  ot  ot 


(4.3) 


where  sa  are  the  dominant  natural  frequencies  or  poles,  va  are  the 
natural  mode  vectors,  ha(s)  are  the  coupling  coefficients,  and  W  is, 
as  it  is  used  here,  the  quasi-entire  function  which  contains  contribu¬ 
tions  from  all  poles  not  included  in  the  dominant  pole  set,  as  well  as 
contributions  from  any  true  entire  function  that  may  exist. 


The  various  quantities  in  Eq.  4.3  can  be  determined  in  the  manner 
14  15 

described  by  Baum.  '  The  poles  are  found  by  a  search  procedure 

(Newton's  method)  as  the  zeroes  of  the  determinant  of  Z  -  Z(s)  .  That 

The  natural 

In  addition,  coupling  vectors,  u^  , 

.  The  coupling  coefficients  are  found 
as  1 


is,  the  poles  satisfy  det(Z  )  -  0  where  Z  -  Z(s  ) 

ot  01  oc 

mode  vectors  satisfy  Z^v  -  0  . 
are  defined  and  satisfy  u^Z^  -  0 


u’v 

h  ( 8  )  - 


U*Z'  V 
a  a  a 


where  V  is  the  value  of  V  at  s  and 

01  ot 


d 

a  3s 


ZL  -  5=  Z(s) 


s-s 


The  coupling  coefficients  are  not  uniquely  defined  and  may  assume  any 

number  of  forms  at  frequencies  away  from  the  poles.  Two  forms  which 

14  15 

will  be  used  for  present  purposes  are  the  "class  1"  form:  ' 
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n\  u^V 

hi1  (s)  -  exp[  ( s  -  s)t  1  J*  — 
*  “  o  uTZ'v 

a  a  a 


where  t  is  the  tum-cn  time  of  the  series  and  the  "class  2"  form: 
o 


T 

u  Z'v 
a  a  a 


The  turn-on  time  is  chosen  as  the  time  at  which  the  impulsive  planewave 
first  touches  the  scatterer. 

Figures  4.1  and  4.2  show  the  geometries  for  the  conducting  sphere- 
capped  cylinder  and  the  conducting  loop  scatterers.  Both  scatterers  are 
bodies  of  revolution  (BGRs)  and  results  were  obtained  with  the  BRFD 
code11  which  uses  a  method-of-moments  solution  of  the  frequency-domain 
electric  field  integral  equation  specialized  for  BORs.  BRFD  can  deter¬ 
mine  the  resonant  poles,  modes,  and  coupling  coefficients,  as  well  as 
the  current  and  scattered  field  responses.  Time-domain  results  are 
obtained  by  inverse  Fourier  transforming  the  frequency-domain  responses. 
The  scatterer  geometry  is  described  by  a  generating  contour  which,  when 
rotated  about  the  z-axis,  sweeps  through  the  surface  of  the  BOR.  The 
surface  coordinates  t  and  4  are  used  to  define  locations  on  the 
BOEl's  surface  where  t  is  the  arc  length  along  the  generating  contour. 
Surface  currents  are  broken  into  t  and  4-directed  components . 

Figure  4.3  and  Tables  4.1  through  ^.4  define  the  dominant  poles 
for  the  two  cylinders  and  two  loops  that  are  used  to  construct  a  pole 
series  which  is  then  subtracted  from  the  total  response  to  find  the 
residual  response  which  corresponds  to  the  quasi-entire  function,  as 
defined  in  Eq.  4.3.  Figures  4.4  through  4.7  show  both  total  and  resid¬ 
ual  t-directed  currents  for  the  L/D  -  10  cylinder  for  6-polarized 


< 
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Table  4.1.  Poles  for  sphere-capped  cylinder  with  L/D  -  10 


s^L/cn 


k 

Real 

Imag 

1 

-0.1404751 

0.822279 

2 

-0.2389305 

1.760744 

3 

-0.3280972 

2.721550 

4 

-0.4203163 

3.68549 

5 

-0.5030614 

4.627963 

6 

-0.5693341 

5.567039 

7 

-0.6114042 

6.515115 

8 

-0.6424545 

7.497416 

9 

-0.6786275 

8.528675 

Table  4.2.  Poles  for  sphere-capped  cylinder  with  L/D  -  30. 


s^L/cn 


k 

Real 

Imag 

1 

-0.1014181 

0.8725323 

2 

-0.1575135 

1.828658 

3 

-0.2035508 

2.803896 

4 

-0.2478629 

3.785218 

5 

-0.2912341 

4.760835 

6 

-0.3279485 

5.724434 

7 

-0.3504116 

6.681224 

8 

-0.3586058 

7.644284 

9 

-0.3629002 

8.620916 
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Table  4.3.  Poles  for  loop  with  b/a  -  10 


skb/c 


k 

Real 

iraag 

1 

-0.2404304 

1.034475 

2 

-0.3971459 

1.989159 

3 

-0.5232771 

2.913365 

4 

-0.6255311 

3.824480 

5 

-0.7108793 

4.730809 

6 

-0.7844273 

5.635828 

7 

-0.8495169 

6.541081 

8 

-0.9083419 

7.447382 

9 

-0.9624029 

8.355284 

Table  4.4. 

Poles  for  loop  with  b/a  -  100. 

k 

skb/c 

Real 

imag 

1 

-0.0958744 

1.045078 

2 

-0.1417687 

2.061150 

3 

-0.1777874 

3.072605 

4 

-0.2085869 

4.081876 

5 

-0.2358296 

5.089798 

6 

-0.2602842 

6.096767 

7 

-0.2823390 

7.103003 

8 

-0.3021917 

8.108644 

9 

-0.3199142 

9.113752 
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incident  plane  wave  from  directions  6  -  90,  60,  30  ,  and  0  degrees  and 
♦  -  0  degrees  at  the  frequency  c/2L  (near  the  first  resonant  fre¬ 
quency  of  the  cylinder).  The  residual  currents  obtained  by  using  both 
class  1  and  class  2  coupling  coefficients  are  shown.  The  incident 
electric  field  magnitude  is  one  volt/meter.  Figures  4.8  and  4.9  show 
the  same  for  frequencies  2  and  3  c/2L  with  incident  angle  6-90 
degrees.  Figure  4.10  shows  ♦-directed  currents  at  the  frequency  c/2L 
and  an  incident  angle  6-0  degrees.  The  phi-di rected  currents,  both 
total  and  residual,  for  other  angles  and  other  incident  directions  are 

almost  identical  in  magnitude  to  the  currents  shown  in  Fig.  4.10. 

4  6 

Figure  4.11  shows  "extended  physical  optics"  currents,  '  so  called  here 
because  they  ignore  shadowing  effects.  That  is,  the  extended  physical 
optics  currents  are  defined  by  J  -  2&  x  H1  over  the  entire  surface  of 
the  scatterer.  To  be  true  physical  optics  currents,  those  currents 
between  6-90  and  180  degrees,  as  well  as  certain  currents  near  the 
ends  of  the  cylinder,  should  be  zero.  While  the  currents  shown  in  Fig. 
4.11  are  for  an  incident  angle  of  6  -  90®  ,  the  magnitudes  of  the  cur¬ 
rents  do  not  change  with  frequency  or  incident  angle. 

Several  features  of  the  currents  are  noteworthy.  First,  the  mag¬ 
nitudes  of  the  class  2  residual  currents  are  almost,  if  not  entirely, 
independent  of  both  the  incident  angle  and  the  frequency.  Second,  the 
residual  currents  show  a  strong  resemblance  to  the  extended  physical 
optics  currents  of  Fig.  4.11  and  a  lesser  resemblance  to  the  true  physi¬ 
cal  optics  currents.  Third,  at  an  incident  angle  of  6-0  degrees, 
the  total  and  residual  currents  are  virtually  identical,  that  is,  the 
dominant  pole  contributions  appear  to  be  minimal  for  this  angle  of  inci¬ 
dence.  The  true  physical  optics  currents  contribute  significantly  more 
to  the  far  fields  than  do  the  currents  of  Fig.  4.11,  since  the  currents 
on  the  front  and  back  sides  of  the  cylinder  tend  to  cancel  for  the  lat¬ 
ter  currents.  However,  in  spite  of  the  resemblance  between  the  extended 
physical  optics  and  residual  currents,  the  contributions  of  the  residual 
and  true  physical  optics  currents  to  the  far  fields  are  comparable  due 
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to  an  imbalance  in  the  residual  currents  between  the  front  and  back 
sides  of  the  cylinder. 

One  should  probably  not  attempt  to  infer  too  much  from  the  resemb¬ 
lance  between  residual  and  physical  optics  responses  noted  above  since 
these  responses  become  decidedly  different  as  the  cylinder  radius 
approaches  zero.  For  the  filament  scatterer,  the  physical  optics 
response,  which  varies  in  direct  proportion  with  the  radius,  is  insigni¬ 
ficant  compared  to  the  actual  response  given  by  Bqs.  2.8  through  2.10, 
which  varies  as  1/log  ka  .  A  refined  local  field  solution  similar  to 
physical  optics  wherein  the  scatterer' s  curvature  is  brought  to  bear  may 
provide  a  more  accurate  description  of  the  residual  currents. 

Figure  4.12  shows  the  total,  residual  and  physical  optics  monosta¬ 
tic  radar  cross  sections  (RCS)  as  a  function  of  frequency  for  the 
cylinder  with  L/D  -  10  .  The  RCS  is  defined  by 

2 

a  -  lim  4nR 
R-*» 


where  E5  is  the  scattered  electric  field  at  the  observation  point,  E1 

is  the  incident  electric  field  at  the  scatterer,  and  R  is  the  distance 

to  the  far  field.  In  Fig.  4.12,  the  RCS  is  expressed  in  decibels  rela- 
2 

tive  to  L  .  The  incident  plane  wave  is  assumed  to  be  linearly  polar¬ 
ized  and  9-polarized.  Figures  4.13  through  4.15  compare  the  6  compo¬ 
nents  of  various  time-domain  electric  far  field  species.  In  each  of 
these  figures,  the  electric  field  is  given  in  units  of  (/?n  R)~^EqL 
where  EQ  is  the  incident  electric  field  magnitude.  The  time-domain 
incident  field  is  given  by 
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cos  0  -  k  sin  6) 

tx  -  p-r/c)2/t 2]  (4.4) 

where  tQ  -  0.14151/c,  is  an  arbitrarily  delay  selected  for  conven¬ 
ience  in  plotting,  p  -  Jc/k,  k  -  |E|,  E  is  the  propagation  vector  of 
the  incident  plane  wave,  and  EQ  is  the  incident  electric  field 

strength.  Note  that  the  incident  field  has  a  Gaussian-pulse  time  his-  < 

tory.  Figures  4.16  through  4.19  show  the  corresponding  results  for  the 
cylinder  with  L/D  -  30. 

From  the  above  results,  the  following  observations  can  be  made. 

When  the  cylinder  is  relatively  fat,  the  residual  RCS  shows  a  strong 
resemblance  to  the  physical  optics  RCS.  As  the  cylinder  becomes  thin¬ 
ner,  the  physical  optics  RCS  becomes  weaker  relative  to  the  residual 
RCS.  This  last  observation  is  consistent  with  the  above-mentioned  fact 
that  the  physical  optics  response  becomes  insignificant  in  comparison  to 
the  actual  response  as  the  cylinder  radius  approaches  zero.  Finally, 
the  class  2  residual  time-domain  responses  closely  resemble  the  aper¬ 
iodic  components  in  the  filament  responses  of  Sec.  2. 

Figures  4.20  through  4.27  show  the  corresponding  results  for  the 
two  loop  scatterers.  For  the  first  loop,  b/a  -  10  and,  for  the 
second,  b/a  -  100  .  These  results  provide  another  example  of  early- 
time  components  in  transient  scattering  responses  and  tend  to  corrobor¬ 
ate  the  above  observations  for  the  cylinder's  responses.  The  excitation 
for  all  the  loop  responses  is  a  ♦-polarized  plane  wave  incident  from 
various  directions.  The  incident  electric  field  is  always  parallel  to 
the  plane  of  the  loop.  Figures  4.21  through  4.23  and  4.25  through  4.27 
show  the  ♦  components  of  various  time-domain  electric  far  field  spe¬ 
cies.  In  each  of  these  figures,  the  electric  field  is  in  units  of 
(.ttiRr^b  .  The  time-domain  incident  field  is  given  by 


,  2LE 

Ei(r,t) - (k 

x  exp[-(t  - 
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(a)  L/t)  -  30  for  6-90  degrees 


(b)  L/D  -  30  for  9  -  60  degrees 

Total  (0) ,  class  1  residual  (A),  class  2  residual  (+),  and 
true  physical  optics  (X)  RCS  for  sphere-capped  cylinder. 


(c)  0-60  degrees 


-2 

0 


10 


(a)  0  -  C  degrees  of  loop  with  b/a  -  10 


(b)  6-30  degrees  of  loop  with  b/a  -  10 

Figure  4.22.  Comparison  of  true  physical  optics  (solid  line)  and  class 
2  residual  (dashed  line)  electric  far-field  time 
histories. 
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(d)  6-90  degree*  of  loop  with  b/a  -  10 
Figure  4.22  {Cant.) 
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(a)  6  -  0  degrees  of  loop  with  b/a  -  10 


(b)  0-30  degrees  of  loop  with  b/a  -  10 

Figure  4.23.  Comparison  of  class  1  residual  (solid  line)  with  class  2 
residual  (dashed  line)  electric  far-field  time  histories. 
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(a)  b/a  -  100  for  0-0  degrees 


2xMfc 


(b)  b/a  -  100  for  0-30  degrees 

Figure  4.24.  Total  (0),  class  1  residual  (A),  class  2  residual  (+),  and 
true  physical  optics  (X)  RCS  for  loop. 
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(c)  b/a  -  100  for  0-60  degrees 


(d)  b/a  -  100  for  0  -90  degrees 
Figure  4.24  (Cont.) 
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(c)  0-60  degrees  of  loop  with  b/a  *  100 


o 


(d)  0  -  90  degrees  of  loop  with  b/a 


Figure  4.26  (Cont.) 


Figure 


(a)  0  ■  Oi  degrees  of  loop  with  b/a  ■  100 


(b)  0  ■  30  degrees  of  loop  with  b/a  -  100 

.27.  Comparison  of  class  1  residual  (solid  line)  and  class  2 
residual  (dashed  line)  electric  far-field  time 
histories. 
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(4.5) 


r 


i 


X  «pl-(t  -  tj  -  p*t/e)Vl|] 


«hm  tQ  -  0.4444Vc  and  the  other  quantities  in  Bq.  4.5  ara  da£inad 
as  in  Bq.  4.4. 
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CONCLUSIONS 


The  results  of  this  effort  indicate  that,  in  a  practical  SEM  re¬ 
presentation  for  general  scatterers,  a  quasi-entire  function  is  required 
in  addition  to  a  pole  series  constructed  from  the  dominant  resonances  of 
the  scatterer.  The  quasi-entire  function  may  be  regarded  as  a  forced 
component  in  the  same  spirit  as  the  familiar  forced  components  in  the 
responses  of  circuits.  The  forced  components  may  find  a  concise 
description  in  a  form  similar  to  the  extended  physical  optics  response 
with,  perhaps,  refinements  to  account  for  the  scatterer's  curvature. 

One  refined  physical  optics  solution  has  already  been  used.*^'*7  The 

2  3 

description  might  also  assume  a  ray-theoretic  form.  ' 

The  authors  cannot,  unfortunately,  apply  the  completeness  argument 
of  Sec.  3  with  absolute  rigor  to  transient  responses  of  general  scat¬ 
terers  since  we  know  very  little  about  the  distribution  of  poles  in  the 
s-plane  for  general  scatterers.  A  possible  rigorous  argument  might  be 

found  by  proceeding  from  the  completeness  of  the  eigen  modes  and  their 

18 

close,  but  complicated  relation  to  the  natural  modes. 

It  is  interesting  to  notice  that  the  forced  components  are  con¬ 
structed  from  poles,  i.e.,  aspect  angle  independent  features.  However, 
resolving  individual  pole  terms  within  the  forced  components  is  well- 
nigh  impossible.  Given  the  nature  of  this  component,  it  seems  doubtful 
that  any  other  aspect  independent  features  can  be  found.  Target  identi¬ 
fication  using  information  within  the  forced  response,  therefore,  would 
almost  certainly  require  knowledge  of  the  target's  orientation  which  can 
frequently  be  obtained  from  radar  tracking  information. 

Realization  of  resonance-based  target  identification  hinges  on  the 
existence  of  dominant  resonances  and  proper  processing  to  deal  with  the 
forced  response.  Ideally,  information  within  the  forced  response  would 
also  be  exploited  for  identification.  This  effort  has  clarified  the 
nature  of  forced  responses  in  transient  scattering  and  has  shown  their 
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dose  relation  to  the  SEN  entire  function.  Anted  with  this  inforeation, 
it  should  now  be  possible  to  explore  procedures  for  exploiting  both 
doednant  resonances  and  the  forced  response  for  target  identification 
and  other  applications. 
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APPENDIX  A 

POLE  SERIES  REPRESENTATIONS  OF  AN  EXPONENTIAL 


Here  it  is  shown  that 


exp(z)  «  PN(z) 


(A.l) 


where 


PN(z)  -  exp(r)  7T  N(Z.  -  Z)' 


(A.  2) 


■rE1’!-"1 

i-1 


(A.  3) 


Taking  the  logarithm  of  Eq.  A.  2  gives 


log  PN  -  y  -  £  log  -i, 


z.  -  z 


i  “[(V1  -  >)2] 

i-l  L 


If  | z |  «  N  and  \z^  -  N|  «  N  for  all  i  ,  then  Eq.  A.4  becomes 


log  PM  «  z 


which  is  equivalent  to  Eq.  A.l.  If/  in  addition, 
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APPENDIX  B 

APPROXIMATION  OF  A  PULSE  FUNCTION 

The  pulse  function  g(t)  -  u(t)  -  u(t  -  1)  ,  where  u(t)  is  the 
unit  step  function,  is  to  be  approximated  by 

N 

f(t)  -  ^  A.t1"1  exp(-Nt) 
i-1 

on  the  interval  (o,®) .  The  approximation  criterion  is  the  integrated 
square  error: 

l2 

g  -  ^  Aiti_1  exp(-Nt)  dt 
i 

Setting  the  partial  derivatives  of  E  with  respect  to  each  Ai  to  zero 
provides  the  system  of  equations, 

N 

j  -  1,  ...,  N  (B.l) 

i-1 

which  can  be  solved  for  the  A^  that  minimize  E  .  In  Eq.  B.l, 

1 

Fj  -  J  t^"1  exp(-Nt)  dt, 
o 
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APPENDIX  C 

METHODS  FOR  CONSTRAINING  PRONY'S  ALGORITHM 


In  the  implementation  of  Prony's  method,  an  nth  order  polynomial 
is  solved  for  its  roots.  The  order  of  the  polynomial,  N  ,  is  the  num¬ 
ber  of  poles  being  sought  in  the  transient  data.  If  the  coefficients  of 
the  polynomial  are  denoted  as  a  then  the  polynomial  can  be  expressed 
as  per  Ref.  19,  as 


“O  +  OjZ1  +  e^Z2  +  ...  +  -  0  (C.l) 

where  is  usually  set  equal  to  unity.  The  N  roots,  Zi  ,  are 
defined  as 


with  s^  being  the  poles  sought,  and  at  being  the  time  step  size  used 
in  the  analysis. 

If  the  value  of  one  or  more  of  the  poles  is  known — that  is,  we 
know  some  of  the  Z ^  — then  the  can  be  substituted  into  Eq.  C.l. 

For  example,  if  s^  or  Z^  is  known,  then  Eq.  C.l  can  be  written  as 

“0  +  “l2!  +  °2^i  +  ...  +  o^zj-  o  (C.3) 

The  N  +  1  polynomial  coefficients  a  are  solved  for  in  Prony's 
method  by  solution  of  the  difference  equation 
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fCHOCOPY  RESOLUTION  TEST  CHART 
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N-l 

—  ~^|| |> i  k  —  0 ,  1#  •••#  Y  “  1(  Y  “  H  ■  N  (C.4) 

P-0 

The  i  and  Ij^  ace  the  samples  of  the  transient  signal  being  ana¬ 

lyzed  and  M  is  the  total  number  of  samples  being  used.  The  value  of 
M  must  be  at  least  equal  to  2N  to  give  N  sets  of  equations  in  the 
N  unknowns  .  However,  if  the  value  of  a  pole  is  known,  then  one  of 
the  N  equations  can  be  Eq.  C.3  and  N-l  equations  of  the  form  of 
Eq.  C.4  can  be  written. 

If  L  poles  are  known  a  priori,  the  L  equations  of  the  form  of 
Eq.  C.3  can  be  written  as 

N-l 

^  ■  — Z^  ,  i  ■  1,  2,  ...,  L  (C.5) 

P-0 

and  r  -  M  -  N  -  L  equations  can  be  written  in  the  form  of  Eq.  c.4  as 
N-l 

^  “^pfK  “  _IK+N  '  k  -  0,  1#  ...»  r  -  1  f  t-m-n-l  (C.6) 

p-0 

Hence,  there  are  still  y  -  M  -  N  total  equations  to  solve  for  the  N 
values  of  ,  however,  the  system  is  constrained  by  the  knowledge  of 
the  location  of  L  poles.  As  is  usually  done  in  Prony's  method,  if 
M  ■  2N  ,  then  the  set  of  equations  is  inverted  and  solved.  If  H  >  2N  , 
then  a  pseudo-inverse  procedure  is  used. 

Using  the  matrix  notation  of  Ref.  19,  if  M  -  2n  and  L  poles 
are  known,  then  we  solve  the  equation 


C-2 


I 


AB  -  C  (C.7) 


where  A  is  a  square  natrix  defined  as 


If  M  >  2N  and  L  poles  are  known,  then  the  solution  takes  on  a 
pseudo-inverse  or  least-squares  form  as 


ATAB  -  ATC  (C.9) 

T 

where  A  is  the  transpose  of  the  matrix  A  and  A  is  now  a  rectangu¬ 
lar  matrix  of  the  form 


